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We study the cyclic dominance of three species in two-dimensional constrained Newman- Watts 
networks with a four-state variant of the rock-paper-scissors game. By limiting the maximal con- 
nection distance R ma x in Newman- Watts networks with the long-rang connection probability p, we 
depict more realistically the stochastic interactions among species within ecosystems. When we fix 
mobility and vary the value of p or Rmax , the Monte Carlo simulations show that the spiral waves 
grow in size, and the system becomes unstable and biodiversity is lost with increasing p or Rmax. 
These results are similar to recent results of Reichenbach et al. [Nature (London) 448, 1046 (2007)], 
in which they increase the mobility only without including long-range interactions. We compared 
extinctions with or without long-range connections and computed spatial correlation functions and 
correlation length. We conclude that long-range connections could improve the mobility of species, 
drastically changing their crossover to extinction and making the system more unstable. 

PACS numbers: 87.23.Cc, 89.75.Fb, 05.50.+q 

The question of how biological diversity is maintained has initiated increasingly more research from multiple angles 
in recent decades p], 0, 0, 0, [H, @] . Mathematical modeling of population dynamics is widely recognized as a useful 
tool in the study of many interesting features of ecological systems. However, the enormous number of interacting 
species found in the Earth's ecosystems is a major challenge for theoretical description. For this reason, researchers 
have built many simplified models to describe the evolution of ecological systems over time d, 0, H, 0, GH [H| • One 
of the simplest cases is of three species that have relationships analogous to the game of rock-paper-scissors (RPS), 
where rock smashes scissors, scissors cut paper, and paper wraps rock. It is a well-studied model of population 
dynamics 0, d Q ES 01, and it can be classified in two ways: three-state or four-state models, depending on 
whether we consider the empty state or not. It is well known that such a cyclic dominance can lead to nontrivial 
spatial patterns as well as coexistence of all three species. 

Recently, Reichenbach and co-workers studied a stochastic spatial variant of the RPS game 0, Qj]] • In their study, 
they run the game with four states: the three original cyclically dominating states and a fourth one that denotes 
empty space. In addition, they introduced a form of mobility to mimic a central feature of real ecosystems: animal 
migration, bacteria run and tumble. They found that mobility has a critical influence on species diversity [ii IBl . 
When mobility exceeds a certain value, biodiversity is jeopardized and lost. In contrast, below this critical threshold 
value, spatial patterns can form and help enable stable species diversity. We shall take this population model as a 
basis to construct a new version of the three-species food chain in the constrained Newman- Watts (NW) networks. In 
the model studied by Reichenbach and co-workers, they consider mobile individuals of three species (referred to as A, 
£?, and C), arranged on a spatial lattice, where each individual can only interact with its nearest neighbors. In this 
study, we introduce some stochastic long-range interactions between elements of the lattice. The stochastic long-range 
interactions occur when there exist long-range connections in the NW networks, mimicking a more real ecosystem: 
e.g., birds can fly, so they c an p rey not only near their nest but also at longer distances from the nest [18j], pathogens 
disperse by air and water pjl [20[, biological invasions related to human influence occur over long distances |2l| , 
and the long-range dispersal of plant seeds is driven by large and migratory animals, ocean currents and human 
transportation [22]. Of course, the long-range interaction cannot be infinite, so we limit the distance of long-range 
interactions to Rmax- That is, the individuals are assigned an interaction distance. For the sake of simplicity, we 
consider that the maximum interaction distance R m ax and the long-range interaction probability p are the same for 
all species. With Monte Carlo (MC) simulations we show that the maximum interaction distance R m ax and the 
long-range interaction probability p play an important role in the coexistence of all three species. 

We consider the four-state RPS model which was described in detail in Refs. 0, [l2|, Q3] • Here, we give a recapitu- 
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FIG. 1: (a) The structure of a constrained NW network. All the long-range connections within the range [l,Rmax] and the 
range of Rmax is 1 < Rmax < L. (b) The characteristic path length as a function of the long-range connection probability p. 
The maximal long-range connection distance Rmax = L/4. 



lation: 

AB A AE, BC A BE, CA A CE. 

AE^AA, BE^BB, CE^CC. (1) 

Here, A, B, and C denote the three species which cyclically dominate each other, and E denotes an available empty 
space. An individual of species A can kill B, with successive production of E. Cyclic dominance occurs as A can kill 
B, B preys on C, and C beats A in turn, closing the cycle. These processes are called 'selection' and occur at a rate 
a. To mimic a finite carrying capacity, each species can reproduce only if an empty space is available, as described 
by the reaction AE — > AA and analogously for B and C. For all species, these reproduction events occur at a rate fi. 
In addition, to mimic the possibility of migration, one can amend the reaction equations with an exchange reaction: 

XY -4 YX. (2) 

where X and Y denote any state (including empty space) and e is the exchange rate. The mobility was defined as 
m = 2eN~ 1 in Ref. [4j], where N denotes the number of sites. From a dynamic viewpoint, the RPS game can be 
described by the mean field rate equations [HI, 03] , 

d t a = a[/i(l — p) — crc], 
d t b = a[/i(l — p) — era], 

d t c = a[/i(l - p) - crb], (3) 

where a, 6, and c are densities of the states A, B, and C, respectively. That is, 

a = N a /N, b = N b /N, c = N c /N, (4) 

where 7V a , A/5, and N c are the number of species of A, B, and C, respectively, p = a + b + c is the total density. These 
equations have a reactive fixed point a = b = c = 3 ^ , which is linearly unstable [l7| • 

The mean field approach does not take into account the spatial structure and assumes the system to be well mixed. 
Therefore, it can only serve as a rough model for dynamic processes. Here, we consider the spatial version of the 
above model in the complex NW network structure [23], and we use the Monte Carlo simulation approach. The two- 
dimension (2D) NW network was constructed as follows: (i) We first built a 2D Lx L (N = L 2 ) regular square lattice. 
So, the total number of connections is 2N. (ii) Then, we randomly chose two sites that have no direct connection. If 
the shortest path length between the two sites was shorter than the maximal distance Rmax, we connected the sites; 
if not, we choose other two sites, until the number of the long-range connections equaled 2pN . Here, the shortest 
path length refers to that we did not take into account the long-range connections. 
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FIG. 2: (Color online) Snapshots of the reactive steady state for m = 4 x 10 -6 , fi = a = 1, and system size L = 1000 (e = 2). 
The long-range connection probability p = 0.08, and the maximal interaction distance increases from Rmax = 2 to 120. 



The above procedure produces a constrained NW network structure as shown in Fig. Ufa). In Fig. [2(b), we plot 
the characteristic path length as a function of the long-range connection probability p. The characteristic path length 
decreases with an increase in the long-range connection probability p. The long-range connection probability p is 
equivalent to the rewiring probability in a Watts- Strogatz network but connections are added without removing any 
of the original ones. So, the modified NW structure is characterized by the probability p and the maximal long-range 
connection distance Rmax- Once the network was built as described above, the evolution of the system over time 
obeyed the following rules (modified from Refs. 0,03 ) : (i) Consider mobile individuals of three species (referred to as 
A, £?, and C), scattered randomly on a square lattice as in Fig. [2(a) with periodic boundary conditions. Every lattice 
site was initially occupied by an individual of species A, £?, or C, or left empty, (ii) At each simulation step, a random 
individual was chosen to interact with a randomly-chosen individual directly connected to it. A process (selection, 
reproduction, or mobility) was chosen randomly with a probability proportional to the rates, and the corresponding 
reaction is executed. 

In the above process, N = L 2 simulation steps constitute one Monte Carlo step (MCS). During one MCS all lattice 
sites had one chance to interact. Over time, the spatial distributions of A, £?, and C species changed from one MCS 
to another, providing the evolution of the system at the microscopic level. 

According to Ref. [l7j], Eq. ([3|) could be cast into the form of the complex Ginzburg-Landau equation (CGLE). 
In accordance with the known behaviors of the CGLE, it was found that the spatial four- state model with diffusion 
leads to the formation of spirals. The spirals' wavelength A is proportional to the square root of mobility 0, [T^ . 
To investigate how the long-range interaction probability p and the maximal distance Rmax affect the behaviors of 
the four-state model, we ran MC simulations of this model in the constrained NW network with periodic boundary 
conditions. All the results that we present were obtained starting from a random initial distribution of individuals 
and vacancies, and each site was in one of the four possible states. The densities of A, B, and C coincided with the 
values of the unstable reactive fixed point of the rate equations ([3|). We considered equal selection and reproduction 
rates, which were set /i = a = 1 [ItJ • So, all four states initially occurred with equal probability 1/4. 

In Fig. El we plot typical snapshots of the reactive steady states for various values of the maximal interaction 
distance Rmax- When Rmax is short, long-range interactions have little effect, and all species coexist. The pattern 
of spiral waves in Fig. [21(a) is similar to the case without long-range interactions. With increasing Rmax-, the spiral 
waves grow in size and eventually disappear for longer enough values of Rmax- When the spiral waves disappear, the 
system becomes a uniform state where only one species exists and the others have died out. This process is similar to 
the result from increased mobility m in the lattice simulation without long-range interactions in Ref. [1] . In addition, 
we computed the extinction probability P ext that two species have gone extinct after 10000 MCS (see Fig. [3]). Fig. [3] 
clearly shows that there exists a critical value R c « 30 in the process of phase transition from coexistence (P ex t tends 
to zero) to extinction (P ex t approaches 1). Of course, the critical value R c depends on the other parameters, such as 
the system size, mobility, and so on. 

To investigate how the long-range connection probability p affects the coexistence, we fixed the maximal long-range 
interaction distance to Rmax — 10 and varied the probability p. The simulation results are shown in Fig. [H and the 
dependence of the extinction probability P ext on p is plotted in Fig. It turns out that p has effects similar to those 
of Rmax on the extinction probability and spiral wave pattern. There also exists a critical value p c w 0.06 in the 
process of phase transition from coexistence (P ex t tends to zero) to extinction (P ex t approaches 1). The critical value 
p c depends on the other parameters in the model, as well. 

In Ref. Q , the authors verified that the spiral wavelength increases with individual mobility and that the wavelength 
is proportional to y/m. They found that there exists a critical mobility M c . When mobility is greater than M c , the 
pattern outgrows the system size, causing loss of biodiversity. 
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FIG. 3: Extinction probability as a function of the maximal long-range interaction distance Rmax- Parameters: L = 200, 
t = 10,000 MCS, \i = a = 1, and m = 1 X 10 -4 . As Rmax increases, the transition from stable coexistence (Pext = 0) to 
extinction (P ex t = 1) sharpens at a critical value R c ~ 30. 



FIG. 4: (Color online) Snapshots of the reactive steady state for m = 1 x 10 4 , ji = a = 1, and system size L = 200 (e = 2). 
The fixed maximal interaction distance Rmax = 10, and the long-range connection probability increases from p = 0.0 to 0.55. 

In this work, we obtain similar results in the case of fixed mobility and variable long-range connection probability p 
or variable maximal interaction distance Rmax- This means that increasing p or Rmax is equivalent to increasing the 
mobility. Although the long-range interaction does not directly change the exchange rate e, it does change the spatial 
structure and leads to faster interactions, particularly for exchange. So, increasing p or Rmax increases mobility m 
indirectly. 

In Fig. [6](a), we plot the dependence of P ext on mobility m in the presence of long-range interactions. With 
increasing mobility m, a sharpened transition emerges at a critical value M c w 1.9 x 10 -4 , which is smaller than 
the value (4.5 ± 0.5) x 10 -4 provided in Ref. [J]. In Fig. [6^b), we also compute P ext without long-range interactions 
(p = 0), holding the other parameters the same in Fig.[6](a). In these conditions, it takes much longer (t = 10N MCS) 
to reach the critical value M c «4x 10 -4 , which approximately coincides with the value (4.5 ± 0.5) x 10 -4 . That is 
to say, in this case the system is more stable than with long-range interactions. 

It is well known that long-range connections change spatial structure dramatically. To learn more information 
about the effect of long-range connections on the emerging spiral patterns, we computed the equal-time correlation 
function qaa (\r — t'\) at r and r' of species A for the system's steady state, which is defined in Ref. [17[, as 



where (...) stands for an average over all histories. 

Fig. [7|(a) plots Qaa obtained from MC simulations. When the separating distance reaches zero, the correlation 
reaches its maximal value. With increasing distance, the correlation decreases and the spatial oscillations appear. This 
oscillation reflects the underlying spiraling spatial structures where the three species alternate in turn. Furthermore, 
the correlation functions could be characterized by their correlation length / corr , which is the length at which the 
correlations decay by a factor 1/e from their maximal value. The value of l corr is obtained by fitting to 
exponentials e~ r l lcorT . This value conveys information on the typical size of the spirals [24| . In Fig. [3(b), we show the 
dependence of l corr on the maximal long-range interaction distance R m ax- The results confirm the scaling relationship 
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FIG. 5: The extinction probability as a function of the long-range connection probability p. Parameters: L = 100, t = 10, 000 
MCS, \i = a = 1, and m = 2 x 10 -5 . As p increases, the transition from stable coexistence (P ex t = 0) to extinction (P ex t = 1) 
sharpens at a critical value p c ~ 0.06. 
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FIG. 6: The extinction probability as a function of mobility: (a) With long-range connections, the transition from stable 
coexistence to extinction sharpens at a critical mobility M c ^ 1.9 x 10 -4 , t = 10, 000MCS\ (b) Without long-range connections, 
M c « 4.0 x 10" 4 . Parameters: fi = a = 1, p = 0.02, R max = L/4. 



l corr oc R m ax for the fixed long-range connection probability p. In addition, it can be observed that the linear fit 
is less good when R ma x around 100. Through extra numerical simulations on the correlation length, we found that 
when Rmax is around or above 100, the system would be at extinction in a high probability. This could affect the 
correlation and correlation length. 

In summary, we studied the influence of random long-range connection on four-state RPS games with NW networks 
based on extensive MC simulations. For a fixed maximal interaction distance Rmax-, as the probability of long- 
range connections p increases, we observe that the spiral waves grow in size and (for larger p) disappear. When the 
spiral waves disappear, the system reaches a uniform state and biodiversity is lost. There exists a critical value p c 
separating coexistence from extinction. Similar behaviors are observed with increasing Rmax fc> r a fixed p. To close 
more ecological realistic model, we also consider the case that XE— > EX occur with a smaller probability where 
X and E (empty place) are not neighboring sites. When p or Rmax increases, we observe that both the size of spiral 
waves grow more slowly and the process of phase transition from coexistence to extinction are slower than before. 

We compared the critical value M c obtained in two cases: with and without long-range connection. It is found that 
M c changes drastically and the systems becomes more unstable if even a weak long-range connection is presented. 
We conclude that long-range interactions could result in improved mobility, and it has dramatic effects on species 
coexistence. This point is also confirmed by the equal-time correlation functions for the system's steady state and by 
the correlation length for different R m ax- 

We are grateful to T. Reichenbach for helpful advice on simulation methods and thank referees' for their construc- 
tive suggestions. This work was supported by the National Natural Science Foundation of China under Grant No. 
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FIG. 7: (Color online) (a) The spatial correlation qaa{t) as a function of r in the reactive steady state, (b) The dependence 
of correlation length l cor r on the maximal long-range interaction distance Rmax. Correlation length is depicted as circle. The 
black line is the linear fit. Parameters: L = 1000, m = 1.2 X 10 -5 , t = 6000MCS, and ji = a = 1 
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